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We describe a calculation of the two-photon decays of heavy vector-mesons, B* -^ B-f^ and D* -^ Djj, 
by using the "heavy meson chiral Lagrangian" . The decay amplitudes are expressed in terms of the 
strong coupling g of the Lagrangian at various powers and the strength of the anomalous magnetic 
dipole fXB'Bj (respectively ^d*d^)- In the charm case we are able to express the branching ratio 
T{D*^ -^ D'^j'j)/r{D*'^) as a function of g only and we expect it to be in the 10^® — 10~^ range, 
depending on the value of g. The determination of fxgtogo^, iJ-B*+B+"f requires a more involved analysis, 
including the consideration of bremsstrahlung radiation for the B*^ -^ 5+77 case. 



1. INTRODUCTION 

The heavy vector mesons D* and B* decay by 
strong and electromagnetic interactions. Since 
the mass difference Mb* —Mb is only 45.78±0.35 
MeV [1], the decay B* -^ Btt cannot occur and 
the main decay of B* is B* -^ Bj. On the other 
hand, the D* — D mass difference barely allows 
the decay to a pion and accordingly the main de- 
cays are D* -^ Dtt, Dj. The strength of the cou- 
plings involved in these processes, go'D-Ki 9d*Dj, 
QB'B-y for the different charge states is not known 
although relative branching ratios for the various 
D* -^ Dtt, D* -^ Dj isotopic spin states are well 
measured [1]. 

In recent years, a theoretical framework has 
been developed for the description of the low en- 
ergy strong interactions [2,3] and of the electro- 
magnetic interactions [3,4] between mesons con- 
taining a heavy quark Q and the pseudoscalar 
Goldstone bosons. The Heavy Quark Chiral 
Lagrangian (HQxL) developed for this frame- 
work [2,3] combines the flavour and spin sym- 
metry of the heavy quark effective theory with 
the SU{3)l (8) SU{3)r chiral symmetry of the 
light sector, allowing the use of chiral perturba- 
tion theory for treating the interaction of heavy 



mesons with low-momentum pions. The symme- 
try of this Lagrangian leads to defining a strong 
coupling g which represents the four couplings 
gB'BTT, gB'B'TT, go'D'-^, go'D-n (as well as cou- 
plings to 77, etc) and a coupling /i which repre- 
sents the anomalous magnetic couplings gB'B'j, 
9B*B-y, 5d*£)*7, go'Dj- This is true in the sym- 
metry limit and corrections will obviously ap- 
pear when one considers deviations from both the 
heavy quark limit and the zero-mass chiral limit. 

At present, there is no experimental determi- 
nation of the strength of any of these couplings. 
On the other hand a variety of theoretical mod- 
els has been advanced for their calculation (see 
[5] for a recent review). As a result, there is a 
rather wide range of theoretical possibilities for 
g and /i, since the predictions of various mod- 
els differ considerably. In the next section we 
shall present a short survey of this situation. One 
might expect a direct determination of go'D-ir 
from the measurement of the decay width of 
D* -^ Dtt] presently, there is only an upper limit 
of rtot(-D*+) < 131keV [6] and the possibility of 
an actual measurement depends on the magni- 
tude of the physical value. It is therefore of ob- 
vious interest to devise independent methods for 
the determination of g and /i couplings. We shall 



describe here our recent suggestion [7,8] of using 
the B* -^ B-f-f and D* -^ D'^-f decays for the 
determination of g and gs'Si- Using the Heavy 
Quark Chiral Lagrangian we shaU express the am- 
phtudes for B* — > Bjj, D* -^ Djj in terms of g 
and II. A detailed analysis will show that the D* 
decay is particularly useful for the determination 
of g while the B*^ -^ 5+77 decay may possibly 
be used to determine the strength of gB*+B+-y 

2. THE THEORETICAL FRAMEWORK 

The decays we are considering, which were 
mentioned firstly in [7], are B*^ -^ 5+77, B*'^ -^ 
B°77, Bf -^ 5O77 and D*+ -^ £"+77, 1^*° -^ 
D°77, Z)*+ -^ Df-f-f. We shaU use for their de- 
scription the HQxL in the leading order of chiral 
perturbation theory. Corrections to it requires 
terms proportional to the mass of the light quarks 
niq, to (I/Mq) and derivatives in the HQx La- 
grangian [9-11]; this question will be commented 
on in the last section. 

The interaction term of HM^L to lowest order 
is [2,3] 

C-' = gTr {Hai^^5A':,H,} (1) 

which defines the basic strong-interaction cou- 
pling g. Ha is a 4 X 4 Dirac matrix, with one 
spinor index for the heavy quark Q and the other 
for the light quark q, 

1 + ^ 



H 



\B:j^-Bj5 , H = -foH^-fo (2) 



and Bl^{v), Ba{v) are the velocity-dependent 
annihilation operators of the meson fields, with 
''^^^aa = 0. a, 6 denote light quark flavours 
(a, 6 — 1,2,3) and the _B*, B fields represent, 
unless otherwise specified, both the beauty and 
charm sectors. A^j, is the axial current con- 
taining an odd number of pion fields, given by 

^M = Wd^c - Ca^C^), C - cM^M/I) with 

M being the 3x3 matrix of the octet of pscu- 
doscalar Nambu-Goldstone bosons. / is the pion 
decay constant, / = 132 MeV. Expanding A^ and 
keeping the first term A^^ = {i/ f)d^A4 we obtain 
the explicit interaction terms from (1) 

23: 
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-Bld^MB^ + H.c. 
f ^ 



■t''. 



(3) 



With the usual definition for the strong B*Bti, 
B*B*Tr [5] vertices 

{7riq)Bivi)\B*{v2e2))=gB*BAq^)q^.4 (4) 



{TT{q)B*{vi,ei)\B*(v2,e2)) 

= .gs'B*7r(g^)ea/3Miye?e29^«r (5) 

the physical couplings are given by the limit 
q'^ -^ m^. Isospin symmetry is expressed by the 
relations 

9B'Btt = 9B' + B»tt+ = —y'^9B' + B+7T« = 

V2gB* + B«7T'> =■ -gB*0B+7r- (6) 

and similarly for the gB'S-n couplings. Compar- 
ing (3) to (4), (5) we obtain these coupHngs ex- 
pressed in terms of 5, and at this point we distin- 
guish between beauty and charm 

gB'BTT = gB'B'TT == {2MB/.f)g , (7.1) 

go'Dn = go'D*. = {2MD/f)g . (7.2) 

Throughout this work, we follow the normaliza- 
tion conventions of [5]. In (7.1), (7.2) mass de- 
generacy for {B*,B) and {D*,D) pairs was as- 
sumed and the use of the pseudoscalar mass value 
in these equations is conventional. 

The incorporation of electromagnetism in 
HM^L is performed by the usual procedure of 
minimal coupling, which leads to the replacement 
of derivatives in the free Lagrangian Cq by covari- 
ant derivatives containing the C/(l) photon field, 
when 

£o = -iTr{Hav''Df,haHb} + 

— Tra^S9^Et . (8) 
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Here S = C^(^) and D^, ^ d^ + V^ with the 
vector current, containing an even number of pi- 
ous, given by V^ = ^(C^S^C + C^mC') 

The interaction vertices of the heavy mesons 
with photons are then [3]. 

(7(fc, e)B{vi)\B{v2)) = eMBivi + V2) ■ e (9.1) 

{j{k,e)B*{vi,ei)\B*{v2,e2)) = 

eMB'{ei-e2){vi+V2)-e . (9.2) 



These vertices are not capable to provide for 
the existing B*Bj interaction and an additional 
gauge invariant term proportional to the electro 
magnetic field F^^ must be added [3,4,11] to the 
Lagrangian Cq + £'"* . It is given as 

£(M) ^ ^Tr{Haa^,F''''Hb5ab) (10) 

which gives when we express Ha from Eq. (2) 

£(A') ^ -e/iF^"^ [iB;+B: + e^^.^rv" {B+ B*^ + H.c.)] 

(11) 
This translates into two "anomalous" vertices for 
the electromagnetic interaction, in addition to 
(91.), (9.2), 

(7(fc,e)5*(«i,ei)|B*(«2,e2)) = 
eiJLMB-{ei-ke-e2-e2-ke-ei) , (12.1) 

{-l{k,e)B{v^)\B*{v2,e2))^ 
-ie^Ms.e^^^^e^rw^e^ . (12.2) 

The new coupling /i is related to the customary 
physical couplings by 

^B*M° = 9b">b»i; Mo'fJ." = 9d'°d»i , (13.1) 



MB'fi^ 



gB*+B+ 



7' 



Md-m^ 



(13.2) 

The values of ^b'Stt, go'Di,, gB'B*n, go'D'Tr, 
gB'B-y, go'D-y, gB'B'j, go'D'-y are not known 
from direct experiments; however, a vast lit- 
erature of theoretical attempts has appeared 
and we expose now succinctly the emerging 
theoretical picture concerning g and /i and 
the respective physical couplings expressed in 
Eqs. (7.1), (7.2), (13.1), (13. 2). It should be 
stressed here that the wide interest in g is due 
also to the fact that its strength appears in the 
expressions for many low-energy electroweak pro- 
cesses, like B -^ Triv, Dg -^ K£jy, B -^ D*t:Iv 
among others (see [5] for comprehensive review). 
One may group roughly the calculations of 
g into four major classes, those based on con- 
stituent quark models, those using relativistic 
quark models, the use of QCD sum rules and the 
effective Lagrangian approach. There are obvi- 
ously different approaches within the same group 
as well as calculations which cannot be classified 
as mentioned. Most calculations obtain, in fact. 



the values of gB'BTr and/or go'DTr- The naive 
quark model result is g = 1 [3] and its modifi- 
cation by the inclusion of chirality [12] brings it 
slightly down to (/ ~ 0.75 — 0.8. The use of rela- 
tivistic quark approaches [13], of which recently a 
new wave of results has appeared gives values of g 
of approximately 0.45-0.65. The QCD sum-rules 
tend to give [14,15] lower values, 5 ~ 0.2 — 0.4. 
We also mention a lattice QCD calculation [16] 
with the result g = 0.42 (4), the analysis of Stew- 
art [11] which uses HQ^L with corrections to or- 
der TUq and 1/Mq and constraints from the D* 
decay branching ratio to obtain g = 0.27^q;q4, 
and a chiral bag model calculation [17] which has 
predicted accurately the observed D* branching 
ratio and gives g = 0.53. We should caution the 
reader on the use in the literature of different def- 
initions for the couplings involved. Thus go'Dir 
defined in [15] is larger by the factor {Mw /Md) 
than the one used, e.g. in Refs. [5,7,12]. More- 
over, in certain works one defines a coupling g (see 
[15]), which is related to g of Eqs. (7.1), (7.2) by 

g = {Md> /Md) (1 + 1^) g, where A is 1/Mq 
correction with a value of 0.7 ±0.1. 

Experimentally [6], there is the upper limit of 
T{D*+ -^ D'^TT+) < 89KeV, which translates into 
g < 0.71. The prediction of the naive quark 
model for this mode is 120-180 KeV, of the rela- 
tivistic quark models 55-80 KeV, of the QCD sum 
rules 10-30 KeV, of Rcf.[17] is 53 ± 3 KeV and of 
Rcf. [11] is 18 KeV. 

The situation concerning /i is very similar. 
Generally, the same models employed for calcu- 
lating g were used also for getting the D* -^ Dj, 
B* -^ B"f modes. Since the relative branching 
ratios for £)*+'0 decays into Dtt vs. £)"'"'°7 are 
known [1], most calculations insist on reproduc- 
ing these ratios; then the calculated absolute par- 
tial widths for D* -^ Dj decays follow the same 
pattern as the calculations of g, i.e. predicting 
rather small width from QCD sum rules [14,18] 
and larger ones from quark models [13], chiral 
approaches [4,17] or potential models [19]. The 
range of variation for the different predictions 
is about one order of magnitude in rate. Par- 
allel predictions are made for the 5*° -^ -B°7, 
_B*+ —f 5+7 decays. To exemplify with typical 
results, the QCD sum-rule approach of Dosch and 
Narrison and of [18] predicts r(B*° ^ B°-f) = 



0.04KeV, T{B*+ -^ B+"f) = O.lKeV, while a chi- 
ral bag model calculation [20] obtains r(i?*° -^ 
B°7) =: 0.28KeV, T{B*+ -^ B+-f) = 0.62KeV 
and even higher values were obtained in certain 
models ([4], [9] and first ref. of [13]). 

3. TWO PHOTON DECAY AMPLI- 
TUDES 

We use now HM^L with electromagnetic in- 
teractions, as detailed in the previous section, to 
calculate the B* -^ B77 process {B stands for 
both beauty and charm). We begin by treat- 
ing the neutral decays only, which are free from 
bremsstrahlung. The calculation is performed to 
leading order in chiral perturbation theory and 
to this order there are no unknown counterterms 
[11,21]. The decay amplitude is given at this or- 
der by [7,8] 
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loop 



(14) 



We shall describe now these various contributions 
to the decay amplitude stressing the couplings, 
without giving the involved detailed expressions 
which can be found in [7]. 

^anomaly rcprcscnts B*° -^ 5"'%"^ B^jj via 
a virtual pion. In the charm case, where the de- 
cay to a physical pion is allowed, we limit the 
s = {ki + /c2)^ variable to be up to 20 MeV away 
from the pion mass. This contribution is propor- 
tional to ag, with a = e^/47r. A^^^^ is given by the 



transition B* 



10 



7 
sertions of /i. Hence A'^j.^^ is proportional to a/ip 



S77, having two in- 

2 



Then we have six classes of diagrams containing 
pions and/or kaons loop, and the photon emission 
may occur from loop or from the external legs via 
the anomalous /i interaction. We shall describe a 
typical diagram for each of these six classes. A-^^J 
describes 5*° -^ {B*+7r^) -^ B^jj with the two 
photons radiated from the virtual charged pion 
and is proportional to ag^. The radiation from 
the virtual B*^ is negligible. A-^^J represents the 
transition B*° ^"B*0"7 -^ (B*+7r-)7 -^ B^^yy 



^S is a simi- 



and is proportional to a/iog^. 
lar class, with the B*^B^y vertex replacing the 
B*^B*^y one in the first step, thus being also 
proportional to afi^g^ . Ay^J^ describes the tran- 



sition S*o ^"(B*7r)"7 ^"S*o "7 -^ SO77 being 
proportional to ag^^Q. Replacing B* by B in the 
loop one gets ^loop which is also proportional to 
ag^^o- Finally, we have i?*° ^■"(-B*"'"7r~)"77 ^■ 
5^77; this describes double radiation by the loop, 
with one photon radiated by the pion and the 
other one due to the anomalous transition fi. This 
amplitude is proportional to afi+g^. 

Calculating the rate of the decay one obtains 
an expression containing 13 terms which depend 
on the products g"iJ,Qfi\. with different powers, a 
and (3 having values between and 4 while 7 has 
values between and 2. 

Turning to the _B*+ -^ ,6+77 amplitude, we 
must consider now also the bremsstrahlung given 
by vertices (9.1), (9.2). Again, we just mention 
here the diagrams involved, explicit expressions 
being given in [22]. Firstly, there are all the dia- 
grams classified as in (14), except that /io and 
/i_i_ are interchanged. However, the major di- 
agrams now are those involving bremsstrahlung 
radiation. Let us denote the additional diagrams 
by i?(^) = Bi'^^ + <^p. Bi'2^ has two classes of 
contributions: (1) the chain i?*+ ^"i3*+"7 -^ 
5+77 where the first vertex has strength e and 
the second e/K+; (ii) the chain i?*+ ^"_B+"7 ^ 
5+77 where the first vertex is e/i+ and the sec- 
ond is e. -BJqqp has similarly a contribution from 
B*+ -^ 7"S*+"^ 7"7r+S*°"^ 77B+ and from 
B*+ ^"7r+B*0"^ 7"B+"^ 77B+. The radi- 
ation from the external legs is due to the (9.1), 
(9.2) vertices. 

In performing the calculations we used for the 
propagator of the vector meson [5] —i{g^'^ — 
v^v'^)/2[{v-k)~A/4:] and for that of the scalar me- 
son i/2[(wfc) +3 A/4], where A ^ Mb'-Mb- We 
also employed physical masses in diagram calcu- 
lations and rates, thus including some features of 
the 1/Mq corrections; moreover, the chiral loops 
we include are themselves of order 1/Mq. 

4. ANALYSIS AND CONCLUSIONS 

The analysis of the various modes involves dif- 
ferent features and must be carried out separately. 
We start with the neutral decay D*'^ -^ D'^jj. 
For this case, the crucial step is to use the ex- 
perimental data on the relative branching ra- 
tios of the strong and radiative I?*+,£)*° decays 



r{D*" ^ D^ttO) : r(D*o -^ D^j) = (61.9 ± 
2.9)% : (38.1 ± 2.9)% and T{D*+ -^ D"7r+) : 
T{D*+ -^ D+7) = (67.7 ±0.5)% : (1.6 ± 0.4)% 
[1]. This allows us to relate (Xq and /i-|_ to g and 
to obtain an expression for T{D*'^ -^ D^ll) in 
terms of g at various powers only [6] 



T{D*^ -^ D°-n) = [2.52xl0-"g2_ 



-5.66xl0""g^ 



^A.76xlO-^q'^+3MxlO-'^"q'^- 



For T{D* 



1.53xlO"V] GcV. 
(15) 



r(D*o ^D%)+r(D*o -^ D 



o^o\ 



we use the same procedure for the relating of /ig 
to g and one obtains r{D*° -^ all) = (2.02 ± 
0.12) X 10~^g^ GeV. This leads to our main re- 
sult concerning the measurement of the value of 
g: we have shown that the measurement of the 
branching ratio 

given by 

Br{D*° -^ D^-n) = 

0.0025 + 0.057g + 4.76^2 + QMg^ + 1.535^ 
2.02 X 105 

is a direct measurement of the strong coupling g. 
In the expressions (15), (16) we assumed relative 
positive signs for 5, /xq, /U+ as indicated by theory 
[11]. However, even if we assume negative rela- 
tive signs for various pairs of the three parame- 
ters (7, /ioj M-i-7 the changes arc very small [6] since 
the main contribution is coming from quadratic 
terms. 

Since g is expected to be in the range 0.25 < 
5 < 1, our result (16) indicates that the branching 
ratio is expected to be 0.17 x 10^'"' < Br{D*^ -^ 
D°jj) < 3.3 X 10"^ The absolute width of the 
two-photon decay obtains from (15) to be 2.2 x 
lO-^eV < r(D*o -^ £'°77) < 6.7 eV for g in the 
0.25-1 range. 

It is of interest to remark here that for 0.1 < 
g < 0.7 the main contribution to the decay rate 
of D*° -^ D'^jj is given by the g^,g^,g'^ terms 
of (15). Out of these terms, the more important 
one is the g"^ one, which accounts for more than 
90% of the contribution for g > 0.25. The 5" 
term gets contributions from the tree term, from 



the A 



0(1) 
loop 



term and from the interference of the 



anomaly with the A\'^ — A^^ terms. However, 
the first one is by far the major contribution, ac- 
counting for about 90% of this term. Hence, the 
tree term emerges as the major contributor to 
the two photon decay D*" -^ D'^j'j for g in the 
"preferred" region 0.2 < g < 0.7. For smaller val- 
ues the anomaly becomes competitive, while for 
g closer to one, the loop terms play this role. The 
differential decay width dT{D*" -^ D°jj)/ds is 
given in [6] for (7 = 0.25 and g — 0.7 and the role 
of the anomaly is visible in the first graph. 

The situation for B*'^ -^ B'^J^, which ampli- 
tude depends also on ^q, /i_|_ and g, is consid- 
erably more difficult to analyze since the help 
of measured relative rates which we had for D* 
decays is unavailable here. The branching ratio 
to the main decay B*'^ -^ B'^^ thus remains de- 
pendent on the three unknown couplings in this 
case. Nevertheless, some help comes from the fact 
that the infiuence of the /i_|_ is very small and 
the process may therefore be analyzed in the re- 
duced parameter space of {/io,ff}- As shown in 
[6], the dependence on the assumed relative sign 
of g/ fiQ is however not negligible now. Moreover, 
the differential decay width does not vary enough 
with the variation of fJ.a,g, in order to provide 
for a reliable determination. It appears that the 
best approach is to wait for the determination 
of g from D* decays. Then the measurement 
of _B*" -^ B°7, both differential and total rate, 
will provide a measurement oi (Iq, — for which 
no other alternative exists. Concerning the ex- 
pected branching ratio, for 0.25 < g < 0.7 and 
40 eV < r(B*0) -^ B°-f) < 1 KeV, Br{T{B*° -^ 
B°-f-f)/T{B*" -^ BO7)} varies between 3.1 x lO"'^ 
and 1.5 x lO"'^ 

Lastly, we refer to the charged decays, £>*+ ^ 
£'+77, D*+ -^ £»+77, B*+ -^ S+77 and let us 
take the latter as an example. Now, the expected 
rate is much larger than in the neutral modes, 
since bremsstrahlung radiation is involved. To get 
an idea of the effect, just taking diagram B^^^l^ 
described in the previous section, one finds a 
branching ratio Br{B*+ —> B'^jj/B*^ -^ B+l) 
of 0.7 X 10^2 foj. ^^^j^^ > lOMeV. This branch- 
ing ratio will be mainly a function of g and /i+ . 
In [22] a detailed analysis is presented of the ef- 
fect of varying the strength of these couplings on 
the branching ratio and differential spectrum of 



i3*+ -^ i3+77. Again, the knowledge of g will 
simplify the problem considerably. 

In concluding, we remark that our calculation 
[6] was performed to the leading order in chi- 
ral perturbation theory and mostly to leading 
order in 1/Mq. Our purpose was to introduce 
the method, to show how one could measure g 
from Z3*° -^ £'°77 decay and to propose a pos- 
sible measurement oi fj,B*-t-B+j, ^b*oboj from the 
analysis of i?* — > i?77 decays. Corrections to 
the leading order should be performed, along the 
methods developed in recent years [5,11,23], but 
these will not affect the qualitative features of the 
method we proposed. 
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